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MULTI-PATH RESISTANT FREQUENCY-HOPPED SPREAD 
SPECTRUM MOBILE LOCATION SYSTEM 



FIELD OF THE INVENTION 
The present invention relates generally to radio communication 
systems and in particular the present invention relates to vehicle or mobile 
locators using frequency hopped spread spectrum communication in a multi- 
10 path environment 

BACKGROUND OF THE INVENTION 
Location of radio transmitters is a technology that has been 
used for many years and in many applications. In general, radio direction 
15 finding is classified into many subclasses using a variety of direction finding 
tools such as antenna loops and antenna anrays. For radio location, the use of 
mobile transmitters and a plurality of fixed base station receivers is generally 
known in the art where the fixed receivers locate the source of the radio 
signal by triangulation of the radio signals. In the alternative, a plurality of 
20 fixed base stations may transmit radio signals and a mobile receiver may 
determine its own position. 

Triangulation is a technique where the fixed base station 
receivers each determine the approximate direction of the radio signal from 
the transmitter. The approximate direction lines are then drawn on a map 
25 where the lines will cross. Since the direction lines are approximate, and 
typically three fixed base station receivers are used, the intersecting area 
between -the three lines is usually shaped like a triangle since the lines rarely 
meet at a single point. The transmitter is then located within the triangle. 
Triangulation of radio signals is fraught with uncertainty. Propagation 
30 interference, signal drop out, insufficient transmitter power, ambiguity in 

bearing measurements received from behind the receiving antennas and multi- 
path reflections are common sources of errors and uncertainty in radio 
direction finding and location. 

Vehicle location systems are known in the art which use a 
35 variety of direction finding and location techniques. One example of a 
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vehicle location system is U.S. Patent No. 4,740,792 to Sagey et al. which 
describes a mobile transmitter and a plurality of satellite receivers. A central 
site determines the location of the mobile transmitters by calculating the 
difference in the arrival time of the spread spectrum transmitter signals at the 
> satellites. This system may be used, for example, for the anti-theft location of 
stolen vehicles or the location individuals carrying the transmitters. 

In the alternative, a plurality of fixed base stations may transmit 
radio signals and a mobile receiver may determine its own position as for 
example in Loran-C or Global Positioning Satellite Systems. One of these 
types of systems is exemplified in U.S. Patent No. 4,665,404 to Christy et al. 
which describes the use of spread spectrum radio transmitters at fixed 
locations and mobile receivers. The receivers locate their own position by 
accurately measuring the time at which the radio signals are received to 
determine the range from the transmitters. 

There is a need in the art for a ground-based vehicle or mobile 
transmitter location system which locates a large number of small, low-power 
transmitters and which operate without confusion or collision. There is a 
further need for a ground-based mobile transmitter location system which is 
immune from multi-path disruption of the radio signals in an urban or 
irregular terrain environment. There is a further need for a ground-based " 
mobile transmitter location system which is capable of carrying messages or 
identification signals simultaneous with the location of the transmitter. 

For the reasons stated above, and for other reasons stated below 
which will become apparent to those skilled in the art upon reading and 
understanding the present specification, there is a need in the art for a 
low-power, spread spectrum vehicle or personal locator system which operates 
effectively in moderate to severe multi-path conditions without collisions of 
simultaneous messages from a plurality of transmitters. 



30 



SUMMARY OF THE INVENTION 
The above mentioned problems with vehicle or personal 
locators and other problems are addressed by the present invention and which 
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will be understood by reading and studying the following specification. The 
present invention is a ground-based vehicle, mobile or personal transmitter 
location system which locates any one of a large number of small, low-power, 
frequency hopping spread spectrum transmitters and which operate without 
5 confusion or collision of the received signals. The present invention is 
immune from multi-path disruption of the radio signals in an urban or 
irregular terrain environment is capable of carrying messages or identification 
signals simultaneous with the location of the transmitter. 

The present invention consists of a central station, a plurality of 
10 base stations and a plurality of mobile transmitters which transmit using a 
frequency-hopped spread-spectrum communication signal. The 
communication signal may be modulated using differential bi-phase shift 
keying, or frequency shift keying, or the like. Each of the plurality of base 
stations include an array of receiving dipole antennas and employs special 
15 algorithms for retrieving very low power frequency-hopped spread spectrum 
signals in a noisy and multi-path environment. The base stations use 
computational algorithms for determining the phase difference between each 
of receiving dipole antennas to determine the direction of the transmitter 
relative to the location of the respective base station. The multiple direction 
20 of arrival angles of the received signal at each base station are corrected based 
on an n-dimensional ambiguity space to locate most probable angles of 
j arrival. The ambiguity space plot is used to eliminate erroneous comparisons 
of dissimilar phases of the incoming signal. The most probable values are 
then tracked over multiple frequency hops and an histogram analysis of the 
25 strongest surviving angles of arrival is performed The two peaks of the 

histogram are used as the two most probable directions of arrival. Each base 
station then communicates the two relative directions of the transmitter to a 
central station where the location of the transmitter is determined by 
tri angulation. The direction of arrival angles from the multiple base stations 
30 are summed in a least mean square approach to find a single direction of 
arrival which is then plotted on a gnomonic projection to correct for the 
curvature of the earth. 
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BRIEF DESCRIPTION OF THH DRAWINGS 
In the drawings, where like numerals refer to like components 
throughout the several views, 

Figure 1 is an overview of the infrastructure of the multi-path 
5 resistant frequency-hopped spread spectrum mobile location system. 

Figure 2 describes the synchronization and message format of 
the outgoing paging signals from the base stations. 

Figure 3 describes the format of the frequency-hopped spread 
spectrum signal transmitted by the remote mobile units. 
10 Figure 4 is a diagram of the receiving antennas of the base 

stations including a plurality of dipole antennas and a single reference omni 
antenna. 

Figure 5 is a three antenna array for an interferometric direction 
finding system 

15 Figure 6 is an ambiguity plane plot of the three antenna array 

of Figure 5. 

Figure 7 is the ambiguity plane for the ratios 5:6 of A/2 the 
three antenna array of Figure 5. 

Figure 8 is an example of multi-path reflections in an urban 

20 environment. 

Fiaure 9 is diagram showing ^angulation of the location of a 
remote mobile unit between three base stations. 

Figure 10 is a diagram of a gnomonic plot of the measured 
bearings of the direction finding algorithm on a sphere. 

25 

. DETAILED DESCRIPTION OF THE PREF ERRED EMBODIMENT 

In the following detailed description of the preferred 
embodiment, reference is made to the accompanying drawings which form a 
part hereof, and in which is shown by way of illustration specific preferred 
30 embodiments in which the inventions may be practiced. These embodiments 
are described in sufficient detail to enable those skilled in the art to practice 
the invention, and it is to be understood that other embodiments may be 
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utilized and that structural, logical and electrical changes may be made 
without departing from the spirit and scope of the present inventions. The 
following detailed description is, therefore, not to be taken in a. limiting sense, 
and the scope of the present inventions is defined only by the appended 
5 claims. 

The present invention is directed to a location system for use in 
locating remote mobile units. The remote mobile units may be located in 
motor vehicles, located on the person of people, carried within containers or 
packages, or any number or variety of mobile carriers. The location system 
10 comprises a plurality of base stations (at least two) which transmit signals to 
the remote mobile units and receive signals from the remote mobile units. 
The base stations may send messages and information to the remote mobile 
units, or the base stations may simply send alert or interrogation commands to 
activate the remote mobile units. Thus the remote mobile units may be 

15 activated locally (for example by the wearer or vehicle) or remotely by the 
base station. The remote mobile units may be used by a person as a reverse 
pager, as an emergency locator or as a communication device. The remote 
mobile units may also be used as part of a vehicle to locate the vehicle if 
stolen, in an accident, for vehicle tracking or as an integral communication 

20 device. 

The remote mobile units receive messages as a standard paging 
device over licensed airwaves using a standard paging infrastructure. The 
remote mobile units, when activated, transmit low-power (less than one watt), 
frequency-hopped, spread-spectrum communication signals. The transmitted 

25 signals from the remote mobile units are received by the base stations which 
are then used to locate the remote mobile units. The base stations are 
. equipped with an array of irregularly-spaced dipole receiving antennas and a 
reference antenna (optional). The direction of the received signals from the 
remote mobile; units is determined by interferometric direction finding using 

30 phase difference comparisons between the various, dipoles of the array 

Unique algorithms using histogram analysis and spatial filtering are performed 
on the received signals at the base stations to eliminate the ambiguity caused 
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by signal interference, signal drop out ambiguity in bearing measurements 
received from behind the receiving antennas and multi-path reflections. 

System Overview 

5 Tbe infrastructure of this acknowledgement paging systems 

serves as the basis of the present invention and the reverse pagers of this 
acknowledgement pager system operate identically in the present system. The 
present invention allows the accurate location of the reverse pagers even in 
noisy and multi-path environments such as those found in urban environments 

10 and irregular terrain environments. 

Figure 1 depicts the major components of the two-way paging 
system in copending U.S. Patent Application Number 08/158,441 entitled 
"LOW-POWER FREQUENCY-HOPPED SPREAD SPECTRUM 
ACKNOWLEDGMENT PAGING SYSTEM" filed November 24, 1993. In 

15 the preferred embodiment of the present invention, all of the components of 
the existing paging system of Figure 1 are used The reverse paging terminal 
1 10 at the central site operates to provide synchronization and messaging 
information through the paging terminal 1 16 to the reverse pagers 100 (also 
known as remote mobile units 100) via direct links to the base stations BSi, ; 

20 BS 2 and BS 3 through ground based radio links (not shown) or through a 
satellite uplink/downlink using a geostationary satellite (not shown). The 
base stations BS b BS 2 and BS 3 include transmit and receive towers 1 13a, 
113b and 113c, respectively and base station terminals 201a 201b and 201c, 
respectively. Terminals 201a-201c are required for producing the accurate 

25 . synchronization information needed to be transmitted to the remote mobile 
units and for local processing of the received messages for. direction finding. 
This synchronization information is used to coordinate the response of 
messaging from the plurality remote mobile units 100 so as to minimize 
collisions within groups of remote mobile units and eliminate collisions . 

30 between groups of remote mobile units. 
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Base to Remote Svnchroniratinn 
Standard paging messages sent from the base stations BS,, BS 2 
and BS3 to the plurality of remote mobile units 100 are, in the preferred 
embodiment, sent as digital data encoded in the POCSAG paging standard 
5 These messages may be used to interrogate the remote mobile units to activate 
the remote mobile unit to allow the base stations to begin the location process. 
Typically the paging channel has a center frequency of 143. 160 MHz, with an 
NRZ FSK data rate of 512 bps or 1200 bps. Other bit rates such as 2400 
baud (bps) are also feasible. Figure 2 describes the POCSAG paging 
10 communications protocol as modified for use by the preferred embodiments of 
the present invention. In the top line.of Figure 2, a greatly compressed time 
line of digital data transmitted according to the POGSAG protocol is shown. 
Batches of messages are transmitted in groups as shown in the details in the 
subsequent lines below the top line of Figure 2. In the second line of Figure 
15 2, a 1.0625 second interval (for 512 baud) is shown in which 544 bits are 
transmitted as a single batch. The batch is preceded by a synchronization 
code word SC as shown in . the third line of Figure 2. This synchronization 
code word is used to get the attention of remote mobile units 100 in the 
geographic locale serviced by the paging terminal. 
20 The synchronization code word within each batch is followed 

by eight frames of digital data Each frame is divided into two portions, an 
address portion and a message portion. The address code word of the 
message of frame 2 of Figure 2 is shown in line 4 while the message code 
word of the second half of frame 2 is shown in-line 5. the address code 
25 word is preceded by. a digital zero followed by 18 address bits, two function 
bits and 10 check bits, The address, code word is followed by an even parity 
bit. The message code word portion of the frame is preceded by a digital one 
followed by 20 message bits which are followed by 10 check bits and a single 
even parity bit. Thus each frame is comprised of 64 bits divided into two 32 
30 bit sections. 

Synchronization of the base station terminal 200 and the remote 
mobile units 100 is necessary to ensure the units 100 are transmitting at the 
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same time that the base, stations are liste ning Synchronization is also 
necessary to coordinate the division of the large number of remote mobile 
units into groups so that members of one group use different frequency 
hopping patterns from members of other groups. Synchronization of the 
5 remote mobile units 100 is accomplished by inserting a special frame into the 
POCSAG data which is used to synchronize the units. 

Synchronization between the reverse paging terminal 110 and 
the remote mobile units 100 is accomplished on two levels. Synchronization 
of the remote mobile units within groups of remote mobile units determine 
10 where along the pseudo random noise code the frequency hops are to be 
followed For example, within a single group of remote mobile units, all of 
the remote mobile units within that group will be synchronized to begin 
transmitting at the same location in the pseudo random noise code list for any 
acknowledgement which may be required Synchronization information is 
15 sent from the reverse paging terminal periodically to the addresses of each of 
the remote mobile units within each group to remind the remote mobile units 
: 100 where along the pseudo random noise code they should be tracking. This 
also enables the dynamic changing of a remote mobile units group 
membership such that if one group is experiencing a large number of 
20 collisions due to simultaneous transmissions, the reverse paging terminal 1 Id 
may re-allocate some of the remote mobile units within that group to new 
groups to minimize collisions. 

Another form of synchronization is required to synchronize the 
remote mobile units to the exact times for transmitting frequencies from 
25 within any of the hops." This fine synchronization information, described 
more fully below, is transmitted as part of the POCSAG codes. 

Referring once again to Figure 2, eight frames of information 
are transmitted in each burst- using the POCSAG format. Remote mobile units 
100 may be assigned to a specific frame within the transmission so that the 
30 remote mobile units, once recognizing the synchronization code word, can 
scan a specific frame for that remote mobile unit's address. Once the address 
is found, the remote mobile unit can detenriine any group changes that may 
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be required to re-allocate that remote mobile unit to a different group. In 
addition, the POCSAG format is used to transmit a fine time synchronization 
code. The fine synchronization code is a transmission of a time pulse at an. 
exact time synchronized to a GPS (Global Positioning System) clock to 
5 synchronize all the remote mobile units 100 for time of transmissioa For 
example, periodically during the day the reverse paging terminal will send a 
synchronization code within the POCSAG code word which is sent at a very 
precise time. In order to ensure that a precise time pulse is sent, the reverse 
paging terminal 1 10 receives accurate time information using a GPS antenna 
10 to receive accurate time of day information. The time used to send the 

synchronization pulse is when the day clock reaches exactly some multiple of 
0.9 seconds in the preferred embodiment In this synchronization information, 
20 bits of information are transmitted to give the accurate time of day 
information. 

15 In each of the remote mobile units 100, the microprocessor 

compares this accurate time pulse which will indicate the exact time of day 
and compare it to its own day clock. The clock within each microprocessor is 
accurate down to a few milliseconds, but the time at which the 
synchronization pulse occurs should have a resolution much finer than that 
20 such as down to 0. 1 milliseconds for time of day. In this fashion, each of the 
microprocessors in each of the remote reverse paging devices can periodically 
realign its day clock to know within a millisecond the exact time. Each 
microprocessor does not actually realign its clock but changes a clock offset 
within memory so that it understands how far off its own internal clock is and 
25 can make the adjustment when using that clock to determine when to start 
transmitting information using the eight frequency hopped spread spectrum 
signal. 

The s\'nchronization pulse is only transmitted every few 
minutes. However, the resolution of the start of the message indicating the . 
30 synchronization pulse is very accurate, it being transmitted at 0.090000 

seconds GPS time after a fixed time of day, such as 12:00 GMT. This GPS 
time is accurate to at least within 100 nanoseconds. 
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An overview of the transmission format of the remote mobile 
unit is shown in Figure 3. The actual transmission of information from the 
remote mobile units 100 is done using Differential Bi-Phase Shift Keying 
(DBPSK) modulation on a frequency hopped carrier of less than one watt. 
5 The transmission of information from the remote mobile units 100 on the 
frequency hopped carrier may also be done using Frequency Shift Keying 
(FSK) modulation. Typically a single transmission consists of 53 hops or 53 
changed frequencies selected from a list of narrow band frequencies. The 
frequency selection is based on a pseudo-random noise code list pointing to 

10 the frequency selection list The synchronization information tells the remote 
mobile unit 100 where along the pseudo random noise code it should be 
synchronized for transmission of its message and the fine synchronization 
information tells exactly the time of day so that the remote mobile unit 100 
knows exactly when to start transmitting the specific frequency so that the 

15 reverse paging terminal 1 10 is looking for that frequency at the same time. 

In operation, 200 frequencies are used by the remote mobile 
unit 100 and the base stations and internally stored in a list numbered Fl 
through F200. For a specific message, 53 frequencies will be used to transmit 
the entire message. These 53 frequencies are selected based on a 1,000 

20 member pseudo-random noise code. 

The use of the accurate synchronization signal periodically 
broadcast via the outbound paging signal enables the remote mobile units to 
use lower accuracy components thus reducing the manufacturing cost of 
remote mobile units. For example, high accuracy crystals to track the, time of 

25 day within the .microprocessor are available with an accuracy of three parts 
per million. Thus, a time drift of approximately three micro seconds per 
second or . 180 microseconds in a minute is the known drift. There are also 
time inaccuracies which are introduced due to the time of transmission 
(variable path length) from the source from the reverse paging terminal when 

30 the synchronization information is sent. By employing crystals which are 

cheaper and have an accuracy of the order 50 parts per million, the amount of 
time-of-day drift normally wouldn't be tolerable. However, by using the 
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synchronization information transmitted on a regular basis from the reverse 
paging terminal, the microprocessor can continually correct its own internal 
day clock so that accurate time of day measurements are always maintained - 
The microprocessor estimates the momentary inaccuracy of the crystal by 
5 tracking the drift across several synchronization transmissions and 

dynamically adjusts for the frequency drift of the crystal and the offset using 
internal offset registers for accurate time of day information. 

Counters are employed within each microprocessor of the 
reverse paging units to compensate for the offset of the frequency based on 
10 the synchronization time information- There are generally two major factors 
which affect the drift in a crystal: temperature and acceleration. Most of the 
drift is due to ternperature, and the remaining drift components are negligible. 
The frequency drift in a crystal due to temperature is very slow, on the order 
of 50 Hz over 10 seconds. During a single day the temperature can change 
15 by 20 or 30 degrees fahrenheit, requiring a time update from the GPS clock 
approximately every five minutes. 

Remote Mobile Unit Transmission Format 
The signal sent from the remote mobile unit 100 to the base 
0 stations is a spread-spectrum, frequency-hopped transmission using differential 
bi-phase shift keying (DBPSK) modulation on the frequency-hopped carrier to 
transmit digital information. The transmission of information from the remote 
mobile units 100 on the frequency hopped carrier may also be done using 
Frequency Shift Keying (FSK) modulatioa The frequency hops are relatively 
> slow, the frequencies transmitted are very narrow and the transmission power 
is extremely small. The maximum peak output power of transmission from 
remote mobile unit 1 00 is limited to less than one watt to allow use of the 
902-928 MHz band in the United States without the need for licensing the 
paging transmitters as allowed by FCC regulations defined in 47 C.F.R 
§15.247, which is hereby incorporated by reference. Those, skilled in the art 
will readily recognize that other- frequency bands and transmissions power 
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levels may be employed depending upon FCC licensing requirements or other 
frequency licensing requirements of other nationalities. 

The use of an accurate crystal to control each frequency of 
transmission is required within each remote mobile unit 100. For example, 
5 high accuracy crystals to transmit the narrow bandwidth frequencies used for 
the frequency hopped transmissions are available with an accuracy of three 
parts per million. At 900 MHz, a 3 ppm drift would place a single frequency 
somewhere within a 2.7 KHz band To tolerate frequency drift due to aging 
and temperature, each individual frequency of the frequency hopped signal is 
10 allocated to a 7.5 KHz band or channel, even though the actual frequency is 
on the order of 200 Hz wide skirt within this 7.5 KHz allocated bandwidth 
Those skilled in the art will readily recognize that by using alternate 
components, the frequency channels (individual frequency of the frequency 
hopped signals) of 7.5 KHz allocated bandwidth may be wider or more 
15 narrow depending upon the overall allocated bandwidth for the system For 
example, 1 KHz or less bands may alternatively be allocated per channel. 

Tests on this invention have shown that by processing the 
received signals at the base stations entirely in the digital domain using the 
combination of unique Fast Fourier Transform algorithms of the present 
20 invention to locate and retrieve the frequency hops and by using a 
combination of unique confidence algorithms with a plurality of error 
correction codes, the receiving base station is able to pull the response 
information from a very low power signal from a distance of up to 45 
kilometers (28 miles) in a flat terrain. In a rather noisy urban environment a 
25 range of 24 kilometers (15 miles) is the norm The information within the 
signals is accurately decoded even though direction of the incoming signals 
may be severely distorted due to multi-path reflections and noise. 

As shown in Table 1, the remote mobile unit message format 
consists of a preamble and the message body spanning a total of 53 frequency 
30 hops. Those skilled in the art will readily recognize that longer messages may 
be transmitted using the preferred embodiment of the present invention, and 
the messages format described here is illustrative but not limiting. Since FCC 
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regulations defined in 47 C.F.R. §16.247 require a minimum of 50 frequency, 
hops, the 53 hop message format described here is considered a short message 
hop. Much longer message hops to transfer more digital data is also 
implemented but not described here. Of course, those skilled in the art will 
5 readily recognize that shorter messages than those described below are equally 
possible for the preferred embodiments of the present invention. The message 
length and number of transmission hops are a matter of design choice. 

The message preamble consists of a predefined code of ones 
and zeros to get the attention of the base unit receiver to begin its FFT (Fast 
10 Fourier Transform) routines to begin pulling the message out of the noise. 
The preamble consists of 165 bits transmitted across 5 hops, that is, 
transmitted using DBPSK (Differential Bi-Phase Shift Keying) or Frequency 
Shift Keying (FSK) on five different frequencies selected from the frequency 
list with the specific frequencies selected based on the PN (Pseudo-random 
15 Noise) Code list stored within the remote mobile unit. The sequence location 
within the PN code that the remote mobile unit will begin to follow is based 
on the synchronized time of day. Within a single hop (a single carrier 
frequency), the carrier phase is modulated 33 times to encode the predefined 
one-zero pattern of the preamble. 
20 The message body follows the preamble and consists of three 

groups of data. Each group consists of 30 actual data bits so that the entire 
message is, in the preferred embodiment of the present invention,; 90 total data 
bits (although other bit length messages may be chosen). The actual (Wa 
encoded within these 90 bits is described above and may be in any convenient 
25 coded format. Those skilled in the art will readily recognize that a wide 

variety of message formats and encoding of the data bits may be used without 
departing from the spirit and scope of the present invention. The encoding 
described here, however, has been proven effective in retrieving the data bits 
buried in background noise with a high degree of accuracy and a low actual 
30 error rate. 



WO 95/14935 'W W PCT/IB94/00371 

• 14 

TABLE 1: Remote Mobile Unit Message Format 

Preamble is 165 bits (33 bits x 5 hops) 
Message is 512 transmitted bits 
(Message is 90 bits actual data) 

£7 = One Frequency Hop . 



10 



I i i r~ 
preamble 
~ i i i i 



i 1 i 1 1 i 1 1 1 1 1 I I i I ' ~ ~~ 1 |. 

Message Body ( Spread Over t 48 Hqps^ ^ j 



Outer Message Coding 
Each of the three groups of message data (30 bits each) are 
15 BCH encoded using a standard 30,63 BCH code and with a single parity bit 
added to form a 64-bit word This encoding decreases the error rate from 10 2 
to 10 -5 . This encoding, documented and understood by those skilled in the art, 
can correct up to 6 errors or detect up to 13 errors. Detection of corruption of 
a data word that cannot be reconstructed will cause the base to request a 
20 second transmission of the acknowledgement message. 

Inner Coding, and Interleaving 
The inner coding of the message will" protect the integrity of the 
message with an error rate as high as 25% Each block of 64 bits of data 

25 (corresponding to a groups of 30 bits and earlier encoded by a standard 30,63 
BCH code) is split into two sunblocks of 32 bits (sub-blocks A and B of 
Table 2), and a reference bit is added to each sub-block to assist the 
differential encoding to provide a reference bit to the DBPSK or FSK 
decoder. The 33 bit sub-blocks are transmitted" over one frequency hop each 

30 and are replicated 8 times so that the 64-bit block traverses 16 frequency 

hops. In transmission, the 33 bit sub-blocks are interleaved to further reduce 
loss of data, as shown in Table 3, where sunblocks A and B of Table 2 
correspond to the first group of 30 bits, sub-blocks C and D, correspond to the 
second group of 30 bits, etc. The total message is 53 hops where each hop is 

35 180 msec in length making the duration of a single message 9.54 seconds. 
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TABLE 2: Interleaving Format for Sub-Mock 

A = 1 reference bit and 32 data bits = 33 bits 
B = 1 reference bit and 32 data bits = 33 bits - 
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TABLE 3: Inner Coding and Interieaving of Sub-Hocks 

A = first 33 bits of 1st block 
B = second 33 bits of 1st block 
C = first 33 bits of 2nd block 
D = second 33 bits of 2nd block 
E = first 33 bits of 3rd block 
F = second 33 bits of 3rd block 
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Those skilled in the art will readily recognize that a wide 
variety of data interleaving may be utilized to effect better error tolerance and 
may be substituted for the interleaving described here. Such alternate 
substitute interleaving means are CIRC (Cross Interleaved Reed Solomon . 
Code) used in CD (Compact Disc) recording media operating either at the 
block level or at the bit level. 



Single Hop Format 
The acknowledgment signals are transmitted by the remote 
35 mobile units 100 in a 1.5 MHz band selected from within the 902-928 MHz 
spectrum. The 1.5 Mhz band is divided into 7.5 KHz channels to provide 200 
channels available in which the frequency hops can occur. Thus, each 
frequency hop is a channel 7.5 KHz wide in which a carrier frequency is 
transmitted. For example, channel one will have a frequency Fl at 
40 902.000000 MHz +/- 3.25 KHz, channel two will have its center earner . : 
frequency at 902.007500, MHz +/- 3 KHz, etc. 
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Each transmit frequency of each hop will thus be centered at 
the approximate mid-point of the assigned channel band; however, due to 
inaccuracies in the remote mobile unit circuits and reference crystals, the 
actual transmit frequencies will vary between units. If high quality crystals 
5 are used to accurately produce the required frequencies, very little drift off the 
center frequency will result In the preferred embodiment of the present 
invention, low cost crystals are purposely employed to keep the per-unit 
manufacturing costs down. This will allow for a lower-cost product sold to 
the user which will increase market penetration Thus, reference crystals are 
10 preferred which have a frequency accuracy of 3 ppm such that at 900 MHz, 
the statistical drift would be approximately 2700 Hz. The crystals center 
frequency within its nominal accuracy also drifts due to aging and temperature 
variations, but this drift is slow compared to the transmission times so the 
drift during a single transmission due to these latter variants is unimportant. 
15 A single frequency hop is shown in Table 4. The 15 

millisecond guard time preceding each hop is primarily a settling time for the 
oscillator circuits of the remote mobile units to allow the internal oscillator 
circuit to lock onto the new frequency between hops. Each hop is transmitted 
at a single frequency in which the phase of the earner is either at 0 degrees 
20 phase or 180 degrees phase in reference to the phase of the reference bit 
immediately following the quiet or guard time. Thus the first bit is a phase 
reference bit followed by 32 data bits exhibiting either zero, phase shift or 180 
degree phase shift to encode the data bits as DBPSK (Differential Bi-Phase 
Shift Keying). In an alternative implementation, each frequency hop may. be 
25 modulated using Frequency Shift Keying (FSK) in which two frequencies are 
used to transmit data bits. One hop frequency may indicate a logical one 
while a second hop frequency may indicate a logical zero. The frequency 
shift is minor and the frequency differential is contained within a single hop 
channel. 

30 Each bit of DBPSK or FSK is a transmission of approximately 

5 milliseconds of the hop carrier frequency either in phase with the reference 
bit transmission or 180 degrees out of phase. Each actual bit is approximately 
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4.7 milliseconds of earner at the hop frequency preceded and followed by 
approximately 0. 15 milliseconds of quiet guard band to reduce discontinuities 
between phase changes which contribute to step- function noise in the 
transmission. 

5 TABLE 4: Single Frequency Hop Format 

Guard time (quiet) = 15 ms 

Single Bit = 5 ms of carrier DBPSK / FSK 

33 Bits plus guard time = 180 ms 

10 



15 



15ms 


5ms 


5ms 


5ms 


5ms 




5ms 


5ms 


Guard 


Ref 


1st 


2nd 


3rd 




31st 


32nd 


Time 


Bit 


Bit 


Bit 


Bit 




. Bit 


Bit 



Frequency Hopping Sequence 
All of the remote mobile units in the market serviced by the 
reverse paging teniiinal for message or location finding use the same pseudo 

20 random noise code to determine the frequency hops. The pseudo random 
noise code is a digital code which is repealing after approximately 1,000 
unique codes, hi the preferred embodiment of the present invention, the 
pseudo random noise code is stored in memory of each of the remote mobile 
units. Those skilled in the art will readily recognize, however, that a linear 

25 feedback shift register could be used to generate the pseudo random noise 
code on a real-time basis instead of using a look-up table which is presently' 
- in the preferred embodiment 

The PN (pseudo-random noise) code list is stored in memory 
and maps to a frequency list. In the preferred embodiment of the present 

30 invention, the PN code list has L000 entries which repeat as a sequence. The 
control means of the reverse paging units continuously maintain a count of the 
proper location within this list for the exact time of day. As described below, 
the time of day for all remote mobile units in the market served by the base 
terminal are periodically synchronized to ensure acknowledgment messages . 

35 are synchronized to transmit the hop frequency at the proper time and to 

synchronize the location within the PN code list that each remote mobile unit 
will use to transmit. 
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The 1,000 member PN code list maps to a 200 member 
frequency list. In order to allow a large number of remote mobile units to 
simultaneously operate in the same geographic market, the remote mobile 
units are divided into groups and the groups are assigned different sequence 
5 segment locations in the same 1,000 member PN list Thus a remote mobile 
unit from group one will begin transmitting a hop at a frequency determined 
from a first location with the PN code, while a remote mobile unit from group 
two may begin transmitting a hop at a frequency determined from a second 
location in the PN code. The remote mobile units from group one and group 
10 two will complete their respective acknowledgement messages in 53 hops. 
Preferably, the sequence of the PN code used to determine the frequencies of 
the 53 hops for the remote mobile unit of the first group will not overlap the 
sequence of the PN code used to determine the frequencies of the 53 hops for 
the remote mobile unit of the second group. More preferable, the frequencies. 

15 chosen based on the non-overlapping segments of the PN code list are 

orthogonal such that the same frequency is never used by two remote mobile 
units belonging to different groups. 

In the preferred implementation, the 1,000 member PN code list 
is divided into 160 hopping sequences. The remote paging units are divided 

20 . into 40 groups with the members of each group synchronized to track the 

same location in the PN code list. The microcontroller of each remote mobile 
unit regardless of its group membership, continuously runs through the 
repeating PN code sequence to stay in synchronization with the base unit and 
all other remote mobile units. Each group of remote mobile units is further 

25 divided into four subgroups such that the remote mobile units within each 

subgroup are assigned one sequence within the PN code list. Although the 53 
hop sequence needed for each acknowledgement transmission may overlap the 
53 hop sequence used by a remote mobile unit in another subgroup, the 
transmission sequences of a remote mobile unit of one group is chosen to not 

30 overlap the 53 hop sequence used by a remote mobile unit in another group. 
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Base Station Design 
As described in the parent patent application, the analysis and 
decoding of the signals received by the base stations from the remote mobile 
units is done almost entirely in the digital domain. The earner frequencies of 
5 the frequency hops are down-converted to a lower frequency in each base 
station and are then digitally sampled The digital samples are then processed 
to locate the phase information of interest for direction finding and message 
decoding. Each base station is constructed with a plurality of digital signal 
processor pipelines which enable simultaneous message decoding and 
10 direction finding of a plurality of simultaneously transmitting remote mobile 
units. 

As shown in Figure 4, the receiving antennas of the base 
stations are constructed using a plurality of dipole antennas 401a, 401b, 
though 401n (generally referred to as 401) arranged in an array, as shown in 
15 Figure 4 and a single reference omni antenna 404. The omni antenna 404 is 
not necessary for the preferred embodiment of the present invention since any 
one of the dipole antennas of Figure 4 may act as a reference antenna for any 
of the other antennas. The omni reference antenna 404 is used as a reference 
antenna because the omni antenna 404 exists as part of the reverse paging 
20 infrastructure upon which the present invention relies and is used for receiving 
the messages from the remote mobile units. 

In the preferred embodiment of the present invention, the 
direction finding (DF) antennas of each base station uses twenty-four sector 
direction finding dipole antennas 401 arranged in four linear arrays as shown 
25 diagrammatically in Figure 4. Each of the four linear airays consists of six 
dipole antennas. The phase between the antennas is measured indirectly - 
each instant of time one of the 24 antennas is selected and its phase is 
compared to the omni antenna 

The switching of the antennas is done every one bit (5 
30 milliseconds) of every frequency hop and is synchronized to the incoming 

data The switching occurs at the 0.7 millisecond gap between the transmission 
of the bits. The antennas 401 are continuously scanned since there could be 
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many remote mobile units transmitting simultaneously from different 
directions. Differentiating between the various simultaneously transmitting 
remote mobile units is capable. If the simultaneously transmitting remote 
mobile units belong to different groups, they will be using orthogonal 
5 frequencies because they are using different sequences in the PN code. If the 
simultaneously transmitting remote mobile units belong to the same group, 
they are nonetheless distinguishable due to the particular unique frequency 
offsets caused by the frequency drift off center frequency of the transmitter 
reference crystal. As described in the above referenced parent application, the 

10 simultaneously received signals may be transmitted within the same 7.5 KHz 
channel during a hopping sequence of remote mobile units from the same 
group, but the base stations are capable of distinguishing and tracking the 
hops of the simultaneously transmitting remote mobile units. This is due to 
the resolution of the Fast Fourier Transforms (200 Hz bins) used in the 

1 5 pipelined receivers, the narrow frequency skirt of the transmitters in the 
remote mobile units and the drift tracking performed in the base stations 
which can distinguish the unique transmission signatures of the simultaneously 
transmitting remote mobile units. 

The direction finding process, described more fully below, 

20 compare the exact phase amplitudes of each of the dipole antennas of the four 

arrays to the phase amplitude at the reference antenna. In the antenna array 
* shown in Figure 4, a circle surrounding the array is divided into eight sectors. 
The actual direction computation is done in the sectors according to Table 5. 

25 

TABLE 5: Sector Allocation on the DF Antenna Array 

30° to 60° - top and left arrays - : ' 
60° to 120° - left array 
120° to 150° - left and bottom array 
150° to 210° - bottom array ' 
210° to 240° - bottom and right array 
240° to 300° - right array . 
300° to 330° - right and top array 
330° to 30° - top" array 



30 



35 
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Those skilled in the art will readily recognize that the whole 
direction finding array does not have to be installed in a single location For 
instance, each linear array of six dipole antennas could be . installed on a . 
different face of roof of a high rise building. In some cases a remote mobile 
5 unit cannot transmit from specific sector. For instance one of the sectors 
faces the sea. In that case the direction finding array could be reduced, and 
only a single direction finding antenna will be installed at that sector. Usually 
the interferometer does not operate in case that a remote mobile unit is very 
close to the base station (less than 1 kilometer). In such a condition, the 
10 other base stations will still perform the direction finding and the central site 
will still be able to perform the to angulation. 

For the reasons exemplified more fully below, the spacing of 
the six dipole antennas in each of the four sides of the array is critical to 
maintaining a low degree of ambiguity in direction finding and multi-path 
15 rejection. The base ratios of the base lines was found by simulations with 
several relative prime numbers and is selected to be 6 : 10 : 4 : 3 : 5 to 
produce a phase gain PG = 56iL The ratio was computed through extensive 
simulations. Although the spacing between the various antenna is critical to 
accurate determination of the incident angle of the carrier, the order of the 
20 spacing of the antennas in the array is not critical. The order chosen here was 
based primarily upon a physical implementation consideration that the array 
; frame be mechanically balanced 

Interferometric Direction Finding 
25 The direction finding technique used in the preferred 

embodiment of the present invention is performed by measuring the electrical 
phase difference between the antenna of the array caused by the angle of the 
incoming wavefront of the spread spectrum frequency hopped signal. In ideal 
conditions where there is no multi-path disruptions, of the wavefront, the 
30 incoming wave from the remote mobile unit is coming in angle of 9 from the 
perpendicular line to the base line on which the antennas are mounted 
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Figure 5 is an example of how the interferometri c direction 
finding technique of the preferred embodiment of the present invention is 
practiced Figure 5 shows, by way of example and not by f limitation, only 
three antennas of the array. The three antennas are used to determine phase 
5 differences of the incoming signals between pairs of the antenna As 

described above, the phase difference may also be measured between each 
antenna and a reference antenna such as the omni antenna 

The interferometer direction finding (DF) technique calculates 
the direction of arrival of the wavefront by measuring the phase of the 

10 incoming wave front in several places in space. This interferometri c method 
is a planar interferometric method which assumes that the antenna array and 
the transmitting remote mobile unit are all in the same relative plane. 
Differences in height of the remote mobile unit and the antenna array are 
ignore and cause litde error in the direction finding of the present invention. 

15 As seen from Figure 5, the three antenna 501, 502 and 503 are 

irregularly spaced such that distance between antennas 501 and 502 is 
great than distance D 12 between antennas 502 and 503. The distance between 
the antennas is selected to be a multiple of a prime number sequence to 
ensure an irregular spacing of all antennas in the array which increases the 

20 probability that an electrical phase difference will always be measured 

between some of the antenna members of the array. The distance between the 
antennas is important to minimize ambiguity errors in measuring the phase of 
the incoming wavefront. If the antennas are spaced too widely apart, there is 
increased ambiguity in detenrdning whether the phase difference of the same 

25 wavelength is being compared or whether the phase of different wavelengths 
are being mistakenly compared. Since the frequency hopped spread spectrum 
signal of the present invention is selected to operate in the 900 MHz band, the 
wavelength X of the incoming signal is approximately 30 centimeters. 

Again, referring to Figure 5, the phase difference between the 

30 two antennas 502 and 503 is know as D I2 and the formula for measuring the 
angle of the incoming signal is : 
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0) {<p 12 +2rcK 12 ) /2ti = D 12 SXN(6) / X 

where - *- 

0 is the angle of the incoming wavefront measured from the 
5 perpendicular of the antenna array baseline, 

<Pi2 is the electrical phase measured between antennas 502 and 
503 

K 12 is an ambiguity factor since the system is capable of 

measuring (p between ± n since the antennas could be 
10 measuring the phase difference of more than one cycle 

X 

D 12 is the distance between the antennas 502 and 503, and 
X is the wavelength of the incoming signal. 

According to this formula, as long as I} 2 is smaller than A/2, the ambiguity 
15 factor K I2 equals 0. Thus the largest antenna spacing in the aixay is selected 

to be smaller that 15 centimeters. 

In order to estimate the inherent instrumental error of the 

direction finding interferometer, the above equation can be written in a 

diflferent form* 



20 



25 



( 2 \ SIN(d) ='((p 12 +27iiq 2 ) / (2tzD 12 /X) 

Assuming that all the errors in estimating 9 is due to inaccuracies in 
estimating cp ? differentiating the above equation produces: * 

( 3 ) . Q (0)=a (cp 12 ) / (2itD 12 COS (6) / A) 



■ where <j is the standard deviatioa 
From this equation, inherent limitations in the system are determined. The 
30 electrical error is divided by the factor 27T D, 2 / X which is termed the phase 
gain (PG) which reduces the error by the ratio between the distance and the 
wave length Also, the error rapidly grows to infinity as 0 approaches 90°. 



WO 95/14935 



PCT/IB94/00371 



24 

Thus the coverage of the present linear interferometer is limited to ± 60° as 
shown by the division of the array into sectors described in above in 
conjunction with Table 5. In order to reduce the error, the phase gain must 
be made as large as possible and each antenna face of the antenna array of the 
5 present linear interferometer array covers only 1/4 of the space. The antenna 
spacing in the array must have a distance smaller than A/2 in order to resolve 
ambiguities. 

Ambiguities in Interferometer Direction Finding 
10 There are two types of ambiguities in interferometer DF: front 

to back ambiguities (the linear interferometer only covers effectively 1/3 of its 
surrounding area) and multiple possible solutions of the interferometer 
equations, in case where the distance between two antennas is greater than 
A/2. This latter ambiguity is caused by measuring the phase difference 

15 between different cycles of the incoming wavefront 

The front to back ambiguity is solved in the present invention 
by using multiple arrays, as shown in Figure 4. In order to cover 360°, the 
present invention uses four sub- arrays of six dipole antennas each, for 
covering 90°. The antennas in each array have a vertical beam width of 120°. 

20 The base station scans all the antennas of the array through the 1:24 antenna 
selector 402 and selects the proper sub-array by comparing the amplitudes 
received from antennas in each array. If two arrays are almost equal in 
amplitude (within 3 dB), then the interferometer pibcess is done on both 
airays and the direction of airival (DOA) of the signal , is computed either 

25 from both or from the proper array, which is selected according to the DOA 
closer to its perpendicular. 

The multiple cycles ambiguity is solved by using more, than 
two antennas in each linear array By way of example, but not by limitation, 
a three antenna array is as shown in Figure 5. The distances between the 

30 antennas that are elected as D 12 =2*A/2 and D23=3*A/2. Using these ratios, the 
equation labeled (2) above becomes: 
(4) SIN{Q) =<p x2 /2n+Ki 2 
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. ( 5 ) SOT(8)=q) 23 /3n+2^ 3 /3 

Figure 6 is an ambiguity plane plot of the three antenna array, 
of Figure 5. Figure 6 shows the phase difference <p 12 plotted against the 
5 phase difference (p^ as an interferometer ambiguity plane plot Figure 6 

shows a family of lines, 601, 602, 603, 604, and 605, parameterized by K 12 & 
K 3 with the same slope created as 9 moves from 0° to ±90°. The following 
formula are used to create the plane plot of Figure 6: 

10 (cp 23 /37t) +2*r 23 /3= (<p 12 /2t0 +;q 2 

The ambiguity resolution algorithm is as follows. The phase 
difference between the antennas is measured as (p 12 and (p^. On the 
interferometer ambiguity plane of Figure 6, find the line represented by K 12 
15 and K s (ambiguity factors) that has the minimal distance, which is the actual 
probability distance, to <p I2 and cp^. After solving the ambiguity, compute 
the exact DOA using the full span of the array i.e. 2.5?. which gives a phase 
gain (PG) of 57L 

One could select ratios of N:N + 1 of A/2, where N is the 
20 number of antennas in the array. This ratio of N:N + 1 of A/2 may be 

selected to be as large as possible in order to improve the phase gain: But as 
N gets larger, the distance between the ambiguity lines gets smaller and the 
algorithm becomes vulnerable to an error in the proper selection of K I2 and 
K^. The result will be a gross error. 
25 If the. measurements distribution of (p is gaussian, then in order 

to compute the gross error probability, the "tail" of the gaussian distribution 
from the distance between the ambiguity lines to "infinity" must be computed 
since the phase is folding after ± 7C . This computation described above is 
only an approximation and is almost true for small values of a((p) which is 
30 true at high signal to noise ratio (SNR) or high signal to interference ratio. 
For white gaussian noise, a ((p) could be approximated by lA^SNR where 
SNR is the signal, to noise ratio. 
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Figure 7 is the ambiguity plane for the ratios 5:6 of A/2. The 
distance between the ambiguity lines became less than hal fcompar ed to the 
ratio of 2:3 which practically increased the gross error to be impractical. 
Thus, the solution is to add more antennas to the anray. Adding more 
5 antennas at different ratios spreads out the ambiguity lines in hyper space 
(rnulti-dimensional space) with dimension of N-l (N - number of antennas in 
the array) and increases the distance between them In the preferred 
embodiment of the present invention, it was concluded that in an urban 
environment six antennas are required in a linear array resulting in a five 
10 dimensional plot of Figure 6 (which is incapable of being shown in a 

drawing). The. following formula are used to create the plane plot of Figure 7 
which is an example of moving the antenna spacing closer together: 



15 



SIN(Q) = (X/D 12 ) * (<p 12 +27t*T 12 ) / (2tz) = (<P 12 /5tO +2K 12 /5 



SIN(e) = (X/D 22 ) * (<p 22 +2*K 23 ) / (2n) = (cp 23 /67t) +K 22 /3 
(cp 23 /67i) +*23/3= (<? 12 /5%) +2K 12 /5 ■ 

Multi-Path Error Sources in Interferometric Direction Finding ~ 
20 The preferred embodiment of the present invention provides 

accurate direction finding and location in urban areas. The major source of 
error is multipath distortion in which the incoming wave front is impaired by 
reflections coming from different directions. Figure 8 is an example of multi- 
path reflections in an urban environment. The unique feature of frequency 
25 hopping enables the interferometer system of the present invention to 

distinguish between the principal reflection and other reflections. Since it is 
not guaranteed that a line of sight exist between the remote mobile unit and 
the base station due to multi-path reflections, the most consistent reflection 
through the hopping message is the one closest to the real bearing. Empirical 
30 resting and observations proved that the bearing error is in practice the size of 
the "urban block" from which the remote mobile unit transmits. In practical 
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tests in an urban environment, the urban block size was about 75 meters in 
average, and up to 200 meters peak error. 

Multi-Path Resistant Dir^ tion Finding Algorithm 
The full direction finding algorithm relates to measurements of 
a linear antennas array having six dipole antennas per side of a four-sided 
array with a single orrmi reference antenna (see Figure 4). The calculations 
are performed on a six antenna subarray after the sub-array is selected to be 
closest to the remote mobile unit as described above in conjunction with 
Table 5. The full interferometer formula for a six antennas sub-array is: 
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where 



<P; is the electrical phase measured between the antennas, 

D; is the distances between the antennas 
X is the wavelength 

0 is the direction of arrival (DOA) 
K; is the ambiguity factor, integer value such that, 

-7t < 27i D ,• * SIN (9) / A. + 2n K ; < 7t ". 
Ej is the additive noise (assumed to be gaussian). 
To solve the equation is (6) above, linear array calculations are 
25 performed to solve the following vector values: 
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X=Q.5 SIN{d) 



The results of these calculations are as follows: 
20 the vector elements (p are limited by: | -0.5< CP; < 0.5 | ; 

(D is the measurements vector and d is known; 
-0.5 < X < 0.5; X is unknown; 

k is unknown; its elements are integer, limited by the elements 

of d 

25 E is the noise vector that is unknown. 

In the calculations, it is assumed that the noise is additive such as: 



30 
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(8) 



. (o 2 i =j 
E{eJ= 0; i e(l,2,3,4,5) 



£'{c i *e i }= i 



2 a 2 i=j 
-a 2 i=j'-l 
-a 2 i=j'+l 
0 otherwise 
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Where : 

<fe is the electrical phase measured at antenna i. 

The solution to this equation is to search among all K e N. the 
20 specific K that minimizes a weighting function in a weighted least square 
procedure (WLS). For every K e N we compute: 

J (9) t i 

25 

The error will be: 
(10) 

30 
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Weighting function: 



no 



dd T R~ r ) T R .J dd T R-i 
sfR^sT 1 ) . I d T R-^ 



10 



Where A is 5 X 5 matrix The rank of A is only 4. In order to 
reduce the amount of real time computation, the following computations are 
performed and the results stored; 

15 O 2 ) A=F T F 

A =AyA s Ay 
[AvA^eigiA) ■ 

f=a£aJ .. ... 

20 Since the rank of A is 4, one of the eigenvalues and thus the 

corresponding eigenvector are 0. Matrices 4X5 and 5X4 are created and 
designated by a ~ over the top of the matrix identifier and thus: 

(B) " . 

25 A=F T F ■ 
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Using these values in the weighting function produces: 



(14) 



<P=F<£ 



K-FK 



10 



The fulT algorithm thus becomes: 



Compute in off-line: 



15 



1 . Compute f . 

2 . For all KEN compute n= fK. 

3 . Compute the matrix E- 



20 



Compute in real time: 



1. Compute: 



25 



30 



2. For all K e N,- compute: 



C{K) 



3. Find the best three (3) Kthat mostly minimize C(K). This 
procedure is done for each hop. 

4. Compute X(K) for those K selected; 
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5. Compute a histogram from all the X(K), using for each DOA 
its distance from the lines (weight) in inverse as weight. 

6. Filter the histogram by passing a boxcar low-pass filter ; £>f 3° 
over the histogram in order to cancel small peaks in the 

5 neighborhood of a real peak. 

7. Find the highest peak and the second to the highest peak. 

8. Report the results. 

Triangulation Process 

10 The input to the location algorithm are reported measured 

directions from several base stations. Each base station reports two directions. 
Each direction is reported with its peak heights which is an indication of its 
validity. The location algorithm projects the direction of the incoming signal - 
on a gnomonic plane tangent to earth, it filters out the false readings, and it 

15 computes the most probable location of the RMU. 

Cassini UTM Projection 
The solution of triangulation on the earth ellipsoid is a highly 
non-linear problem. In order to simplify the problem, we transfer the 
20 observed data to a plane called the gnomonic plane as shown in Figure 10. 

The geographic location of the base station on earth is transferred to an X-Y 

v ■ 

* position on the gnomonic plane. Each measured bearing is also transferred to 
the gnomonic plane with a small adjustment. The estimation process is done 
on the plane. The estimated location of the remote mobile unit is not 

25 transferred back to an earth ellipsoid since the same plane is used for graphic 
.display of maps using a Microsoft Windows program or the like. 

The Cassini UTM projection is widely used for mapping. Its 
advantages are that the mapping are in the Cartesian coordinate system such 
that the E-W and N-S lines are orthogonal. The projection maintains 

30 approximately both the directions related to the north and the distances at 
radius less than 50 M. The projection is very convenient for graphic display 
such as in the Microsoft Windows format. 
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The main distortion of the UTM projection is caused by the 
convergence of the meridian toward the north- There are two corrections that 
have to be made: .-ar ?• ^ 

5 1. Correction of the map north to the geographical north, which is 

a bias correction; and 
2. Small correction of the DOA (direction of arrival of the signal) 
according to the location of the base station. 

10 In actual observations, the correction was actually measured 

together with the calibration and measurement of the pointing of the direction 
finding array toward the north. The calibration was performed relative to 
known reference point locations in the covered area To improve system 
accuracy, a reference transmitter having a known exact location may be used 

15 to correct the system parameters. , 

Direction Finding Fixing on a Plane 
The fixing process input are the two DOA received from all the 
base stations. There are several tests performed on the results received from 
20 the BS: 

If the peaks are too low, both measurements are ignored. 

If the ratio between peaks is larger than 4, the second peak is 

ignored. 

The CS (Central System) computes the fix by using all the combinations from 
25 the DOA received from the base stations. Note that within one message there 
could be several DOAs. 

Hie fixing process is an LMS (least mean square) process in 
which the most probably point is found The most probable point is defined 
as the point on the plane that its sum of the squared distances from all the 
30 DOAs is minimal. The process is as follows: 

The distance of a point (X,Y) from a DOA is: 
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dtiXrY) =(y-r i )cos(6 i ) - (X-X ± ) siniOJ 



■where Y f ) are the coordinates of the^ase stationonumber i. 
d -, (X , Y) - distance of (X , Y) from the DOA number i. 
5 & - The angle of DOA number i from the north 

(X , Y) - a point on the plane 

(15) could be written as: 

10 (16) j^sinfO^ cosfBJ =X sin (6^ -r cos(8 i ) +d i 

Assuming that M inputs from the base station, (i € {1, ,M}), 

equation (16) is written in a matrix form 



15 (17) 



X d sinidt) -Yv cos Oj) 



js^ sin(6„) -Y M cos(6 w ) 



sin(9 i ) -cos (e_£) 



sin(0 w ) -cos (6^) 



AT' 



r ■ ■ 



20 



where M - number of bases, as shown in Figure 9. 
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(18) 



10 



15 



'X ± sin^) -Y i cos (6J 



X M sin(0„) -Y M cos (6 M ) 
sin(6 i ) -cos(6 i ) 



sin'O^) -cos(d M ) 



D= 



Z=AC+D 



Where: 

z is a vector of M elements. 
20 A is a matrix of M rows, 2 columns. 

c is a vector of M elements. 
D is a vector of M elements. . .* 
The purpose is to find z such that D T D will be minixriized. The 
algorithm is LMS, as described above, and is performed by differentiating 
25 D T D relative to c and comparing the result to 0. 
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Complete Fixing Algorithm 
Compute all intersections between all probable DO As using 
(19). f 

Perform two dimensional filtering over all intersections, taking 
into account the power of each intersection according to the 
DO As power that created that intersection. The filter size is 
900 X 900 meters. 

Find the an intersection from the DOA that are in the filter. 

Operation S ummaiy 
In summary, and in referring once again to Figure 6, the slope 
of the lines is 2/3 and the algorithm finds the closest line to the observed 
point 600. The most probable actual point lies on the line 601. Figure 7 
shows the use of more ambiguity fines to get better resolution of the observed 
15 point 701. The actual probable location of the observed point would be on 
the line at position 702 where the slope of each of the ambiguity lines is the 
distance between the antennas. Thus a closer spacing of the antennas shown 
at Figure 7 produces less ambiguity as to the observed versus actual angle of 
the DOA 

20 The ambiguity lines shown in Figures 6 and 7 are for three 

antennas. In reality according to the matrices described above, six antennas 
are used in a five dimensional space which is unable to be drawn in a figure. 
In a five dimensional space, the observed point is located, and the nearest 
lines in the five dimensional space are located for the most probable actual 
25 DOA for each based statioa By spreading out the spacing of the antennas, 
more ambiguity lines as shown in Figure 7 are used such that antenna spacing 
is best spread out to minimize ambiguity. 

In the based station, the observed point 701 representing DOA 
0 is compared to the closest three lines in a five dimensional space. A 
30 histogram is then generated based on the three lines closest to the observed 
point. Then a gaussian distribution of the observed 0s is computed and only 
the actual points 702 and five dimensional space are picked which lie within a 



5 



3. 



10 
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single standard deviatioa . A histogram typically produces two peaks for most 
probable observed values. Histograms that produce more than two peaks are 
filtered so that only the two strongest peaks are taken as potential values? 
These two most probable actual values of the DOA 0 are then sent to the 
5 central site for triangulation to determine the position of the transmitter. 

Those skilled in the art will readily recognize that the benefit of 
the histogram analysis may be achieved through a less computation intensive 
operation known as K/m where the DOA 0 most often occurring in the 
analysis is the most probable DOA 0. A complete histogram need not be 
10 computed if instead a statistical threshold is applied in which the number of 
occurrences of 0 exceeds a threshold value. As soon as the threshold value 
is reached, further analysis is halted since the most probable 0 has been 
found. 

At the central site, the two most probable 0 for the two 
1 5 histogram peaks around a plus or minus three degree range are recorded for 
two, three or more base stations. Thus 2 1 most probable points are received 
by the base station whereas n is equal to the number of stations. Thus if four 
stations are reporting, two most probable DOA 0, J6 observed points are 
received at the central starioa A spacial filter is used with a 900 by 900 
20 meter space to determine which of the points cluster within the most probable 
location. Within this 900 square meter space, the most probable location is 
selected Additional computation can be made by overlaying the 900 by 900 
square meter space onto a map of an urban area and, for vehicle location, 
determining that the vehicle could only be on a street and not within a 
25 building. In actual observations, the actual location of the remote mobile unit 
is within a 77 square meter area 

Thus the multipath problem is solved through the use of a 
multidimensional ambiguity space to eliminate observed DOA values which 
are least probable. The use of a plurality of arrays will produce a number of 
30 observed values within the n-dimensional ambiguity space, some of which are 
due to the actual direction of arrival and some of which are due to multipath' 
reflections. Since the observed values in the preferred embodiment of the 
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present invention are the result of multiple frequencies in a frequency hopped 
spread spectrum system, observed values and most probable values will be 
calculated at various frequencies. Different frequencies observed qyer~a 
period of time may be more or less susceptible to the multipath problem In 
5 addition, for a moving remote mobile unit, the doppler shift will be observed 
but minimized over time since only the most probable values due to the 
strongest signals will survive a histogram analysis. Thus the histogram 
analysis should be performed over a number of frequency hops over a longer 
period of time until the multipath signals statistically drop out. Thus only the 
10 true direction signals will survive the histogram analysis and will be passed to 
the central station for final analysis. 

CONCLUSION 
Although specific embodiments have been illustrated and 
15 described herein, it will be appreciated by those of ordinary skill in the art 
that any arrangement which is calculated to achieve the same purpose may be 
substituted for the specific embodiment shown. This application is intended 
to cover any adaptations or variations of the present invention. Therefore, it 
is manifestly intended that this invention be limited only by the claims and 
20 the equivalents thereof 
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WHAT IS CLAIMED IS: 

1- A mobile transmitter location system, comprising: r - ~ 

a mobile transmitter having a unique identifier, and having a 
5 transmitter for sending a frequency-hopped spread spectrum carrier signal 
containing the unique identifier, 

a first base station having an array of inregularly spaced dipole 
antennas operable for receiving the frequency-hopped spread spectrum carrier 
signal and capable of determining a first direction of origin of the frequency- 
10 hopped spread spectrum carrier signal; 

a second base station also having an array of irregularly spaced 
dipole antennas operable for receiving the frequency-hopped spread spectrum 
carrier signal and capable of determining a second direction of origin of the 
frequency-hopped spread spectrum carrier signal; and 
15 a central station connected to the first base station and the 

second base station, and operable for. 

1.) receiving the first direction of origin of the 
frequency-hopped spread spectrum carrier signal from the first 
base station; 

20 2.) receiving the second direction of origin of the 

frequency-hopped spread spectrum carrier signal from the 
second base station; 

3.) triangulating the first direction with the second 
direction to determine the location 

25 

2. The mobile location system according to claim 1 wherein a reference 
station operates periodically to transmit a reference frequency-hopped spread 
spectrum carrier signal of known origin and location to calibrate the first and 
second base stations and the central station 

30 

3. The mobile location system according to claim 1 wherein the first base 
station and the second base station both include a reference antenna such that 
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the phase difference of the frequency-hopped spread spectrum earner signal 
between each dipole of the anray of irregularly spaced dipole antennas and the 
reference antenna is determined. *r ?- *z 

5 4. The mobile location system according to claim 1 wherein the first base 
station and the second base station both operate to reduce the ambiguity of the 
direction of arrival of the incoming signal by calculating a multidimensional 
ambiguity space having lines indicating the only possible angles of the 
direction of arrival of the signal and detennining the probable points of these 
10 lines which are closest to the observed value. 

5. The mobile location system according to claim 4 wherein the first base 
station and the second base station both generate a histogram of the probable 
values and sending the two most probable values to the central station 

15 

6. The mobile location system according to claim 1 wherein the central 
site received the most probable values from the first and second base stations 
and applies a spatial filter to only use the most probable values which lie 
within a preselected distance of the average location of all most probable 

20 . values. 

7. A method of locating a mobile transmitter location^ comprising 
the steps of: 

receiving at a first base station a radio frequency signal on a 
25 plurality of irregularly spaced antennas; 

receiving at a second base station a radio frequency signal on a 
plurality of irregularly spaced antennas; 

calculating the observed angle of arrival of the radio frequency 
signal by comparing the phase difference of the radio frequency signal 
30 between pairs of the plurality of irregularly spaced antennas; 
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determining the most probable angle of arrival of the radio 
frequency signal by comparing the observed angle of arrival to ambiguity 
lines in an n-dimensional space where the number of antennas isf irH; 

calculating an histogram of the most probable values and 
5 selecting the strongest most probable values; 

applying a spatial filter to the strongest most probable values 
from each base station and selecting the location based on the closest 
grouping of the strongest most probable values. 
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